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Our understanding of various states of matter usually 
relies on the assumption of thermodynamic equilibrium. 
However, the transitions between different phases of mat- 
ter can be strongly affected by non-equilibrium phenom- 
ena. Here we demonstrate and explain an example of non- 
equilibrium stalling of a continuous, second-order phase 
transition. We create a superheated atomic Bose gas, in 
which a Bose-Einstein condensate (BEC) persists above 
the equilibrium critical temperature, T c , if its coupling to 
the surrounding thermal bath is reduced by tuning inter- 
atomic interactions. For vanishing interactions the BEC 
persists in the superheated regime for a minute. However, 
if strong interactions are suddenly turned on, it rapidly 
"boils" away. Our observations can be understood within 
a two-fluid picture, treating the condensed and thermal 
components of the gas as separate equilibrium systems 
with a tuneable inter-component coupling. We experimen- 
tally reconstruct a non-equilibrium phase diagram of our 
gas, and theoretically reproduce its main features. 

Non-equilibrium many-body states can persist for a very 
long time if, for example, a system is integrable, the transition 
to the lower free-energy state is inhibited by an energy barrier, 
or the target equilibrium state is continuously evolving due to 
dissipation. Ultracold atomic gases offer excellent possibil- 
ities for fundamental studies of non-equilibrium phenomena 
ifT UTZl and have been used to create counter-intuitive states 
such as repulsively bound atom pairs [3| and Mott insulators 
with attractive inter-particle interactions [ 10 1 . 

Our superheated Bose gas is reminiscent of superheated 
distilled water, which remains liquid above 100 °C. Specif- 
ically, as the temperature characterising the average energy 
per particle and the populations of the excited states rises 
above T c , the cloud remains in the partially condensed phase, 
which in true equilibrium should exist only below T c . How- 
ever, there are also important differences. Boiling of water is 
a first-order phase transition and is stalled in clean samples by 
the absence of nucleation centres. In that case the transition is 
inhibited by an energy barrier. For a second-order phase tran- 
sition such a barrier does not exist and the superheating we 
observe is a purely dynamical non-equilibrium effect, which 
arises because different properties of the system evolve at dif- 
ferent rates. In this respect our gas also bears resemblance 
to the long-lived non-equilibrium spin structures observed in 
spinor condensates J4][7), pre-thermalised states in quenched 
one-dimensional (ID) Bose gases 11211 and supercritical su- 
perfluids predicted to occur in quenched 2D gases [ 13 1. In all 
those cases, however, non-equilibrium states are observed due 
to the system's slow approach to a true equilibrium. Here, the 
system actually evolves away from equilibrium. 



In Fig. [T] we summarise the basic idea of our experiments 
and the key concepts needed to understand them. In an equi- 
librium gas, a BEC is present only if T < T c , where T c de- 
pends on the total particle number N, or equivalently if the 
chemical potential is /i > /i c . In a standard experiment, after 
a BEC is produced, it gradually decays because T rises, due 
to technical heating, and/or T c decreases, because N decays 
through various inelastic processes. As T/T c increases, elas- 
tic collisions redistribute the atoms between the thermal and 
condensed components, aiming to ensure the equilibrium par- 
ticle distribution. The BEC atom number, A*o, can therefore 
decay in two ways: (1) by direct inelastic loss, and (2) through 
elastic transfer of atoms into the thermal component. Here we 
reduce the rate of the elastic particle transfer by tuning the 
strength of inter-particle interactions, characterised by the s- 
wave scattering length a. This protects the BEC deep into the 
superheated regime, where No > even though T > T c . 

We can understand our observations within the two-fluid 
picture outlined in Fig. |TJb). Here we treat the thermal and 
condensed components as two coupled sub-systems with atom 
numbers N' and Nq, chemical potentials p! and fj,Q, and in- 
stantaneous per-particle inelastic decay rates T' and To, re- 
spectively. In equilibrium p! = p ; note that p is defined 
only if No > 0, so po > p c . 

The two components are coupled in two ways, both depen- 
dent on the scattering length a. First, the local "kinetic" ther- 
mal equilibrium between the collective excitations in the BEC 
(phonons) and the thermal bath is ensured by Landau damp- 
ing, the rate of which is oc y/a lfT4T - fl6l . Second, the global 
"phase" equilibrium (i.e., the equilibrium condensed fraction 
Nq/N) is ensured by the elastic scattering with a rate oc a 2 . 
Crucially, due to the different scalings with a, we find a large 
parameter space where the two components can be consid- 
ered to be in local kinetic equilibrium while the system is not 
in global phase equilibrium. In other words, the two compo- 
nents are at the same temperature, but have different chemical 
potentials. 

In our optically trapped 39 K gas ifTTl . we control a by an 
external magnetic field tuned close to a Feshbach resonance at 
402 G 1 18 1, the dominant source of V and r is spontaneous 
scattering of photons from the trapping laser beams, and Tq 
has an additional contribution from three-body recombination. 

The key steps in our experimental sequence are summarised 
in Fig. |TJ C )- We start by preparing a partially condensed gas 
in the \F,mp) = |1, 1) hyperfine ground state by evaporative 
cooling at a = 135 do, where ao is the Bohr radius IfTTl . We 
then reduce a (over 50 ms) and follow the subsequent evolu- 
tion of the cloud, probing the atomic momentum distribution 
by absorption imaging in time-of-flight expansion. Reducing 
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FIG. 1 : Creating and understanding a superheated Bose-condensed gas. (a) In equilibrium, a BEC is present if T < T c or equivalently /i > fi c 
(here /? = 1/(&bT)). The arrow indicates the cooling trajectory along which a BEC is produced. The insets show measured momentum 
distributions, with the condensed component indicated in red. (b) Two-component picture. The thermal and condensed components have 
chemical potentials fi' and no, and inelastic decay rates T' and To, respectively. The net flow of particles between the two components, k, 
depends on y! , fj-o and the scattering rate oc a 2 . In equilibrium fx' = fj,o = fi cq . The Landau damping of the collective modes in the BEC has 
a rate oc y/a. (c) Time-sequence of the experiment. Reducing a after preparing a BEC reduces the coupling between the two components and 
extends the condensate lifetime. In the superheated regime // < /i cq < fi c but fio > He- 



ll (at constant No) initially reduces /xrj below p! ifTTI . but sub- 
sequently fi decays slower. 

In Fig. [2] we quantitatively contrast the equilibrium evolu- 
tion of a cloud at a = 83 ao and the non-equilibrium evolution 
at 5 ao. In both cases we start at time t = [see Fig.[TJc)] with 
N « 2 x 10 4 and N w 2 x 10 5 at T 160 nK. In both cases 
T c decreases at a similar rate due to similar TV-decay. At 5 ao, 
the temperature rises faster due to less effective evaporative 
cooling at a fixed optical trap depth. 

Whether the gas is in equilibrium or not, it can always 
be characterised by two extensive variables, the total particle 
number N and energy E. We measure these quantities by di- 
rect summation of the momentum distribution and its second 
moment. To measure N we count the atoms within the cen- 
tral peak rising above the smooth thermal distribution. 

From the measured N(t) alone we calculate the equilibrium 
T c (t) ED- From N(t) and E(t) we calculate the equilibrium 
intensive thermodynamic variables /x cq (i) and T oq (i), and the 
equilibrium number of condensed atoms, N^ft) |[T9l - l2"Tl ; in 
these calculations A^ q > if and only if T cq < T c . For com- 
parison, we also directly fit a temperature T* to the wings of 
the momentum distribution. Additionally, supposing only that 
the two components are separately in equilibrium, from the 
measured Nq and N' we calculate [iq and //. (For theoretical 
details see Supplementary Information 1 16 1.) 

At 83 a [Fig. |2|a)] we find excellent agreement between 
the measured No and the iVQ q predicted without any free pa- 
rameters. The BEC vanishes exactly at the equilibrium "crit- 
ical time" t c (dotted green line), at which T cq = T c . Note 
that the dashed red lines show the experimental bounds on the 
time t when the BEC vanishes. The separately calculated fio, 
fi' and fi cq are all consistent and we have also checked that the 
fitted coincides with the calculated T oq . All this gives us 
full confidence in our equilibrium calculations. 

At 5ao [Fig. |2jb)] we observe strikingly different be- 
haviour. The BEC now survives much longer than it would 
in true equilibrium; t — t c « 40 s. We also see that /xq and \J 



diverge from each other for t > t c , so the system is moving 
away from the global phase equilibrium rather than towards 
it. The observed superheating can thus not be understood as 
just a transient effect. (Note that /io — is always very small 
due to weak interactions.) 

At 5 ao the gas is not in global phase equilibrium, but it 
is still a good approximation to view its two components as 
two equilibrium sub-systems at a same (kinetic) temperature, 
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FIG. 2: Equilibrium vs. non-equilibrium BEC decay, (a) At a = 
83 ao the cloud is always in quasi-static equilibrium. The measured 
iVo is in excellent agreement with the predicted A^Q q and vanishes 
when T cq = T c \ the three separately calculated chemical potentials, 
Ho, fi' and n cq , all agree with each other. The dotted green line marks 
the equilibrium critical time, t c , and the dashed red lines show the 
experimental bounds on the time i when the BEC actually vanishes, 
(b) At 5 ao, the BEC persists in the superheated regime (T oq > T c ) 
for t — t c w 40 s. 
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FIG. 3: Quenching the superheated Bose-condensed gas. Solid sym- 
bols show the evolution of No at a = 3 ao, the green solid line shows 
and orange shading indicates the superheated regime. Open 
symbols show the rapid decay of the BEC after it is strongly coupled 
to the thermal bath by an interaction quench to a = 62 ao at time 
t q . We show two experimental series in which t q = 20 s (black) and 
30 s (red). 

as in Fig. |TJb). We have checked that the momentum dis- 
tribution in the non-condensed component is still fitted well 
by a thermal distribution, with T* always within 10 % of the 
calculated T eq (N,E) (see Supplementary Information ifTrjlO . 
For the BEC, in a weakly interacting gas the equilibrium rela- 
tion no (No) relies on the macroscopic occupation of a single 
quantum state jl22ll . rather than on global equilibrium. More- 
over, even for the lowest-energy collective modes we estimate 
the Landau damping time to be < 1 s |[T4l - tT6l . i.e. much 
shorter than the characteristic time scale of our experiments. 
Thus, while this distribution is not directly measurable, we 
expect the distribution of collective excitations in the BEC 
to be characterised by a temperature To that is also close to 
T cq s» TA 

These conclusions hold for any a > 1 ao 11161 . Exactly at 
a = our theoretical picture does break down, since the Lan- 
dau damping rate vanishes and the BEC has no equilibrium 
features; the two components are simply completely decou- 
pled. Bearing this small caveat in mind, from here on we refer 
to T ~ T' ss T cq simply as the temperature of the system 
T. 

If superheated water is perturbed, e.g., by sprinkling some 
salt into it, it rapidly boils away. Here, an analogous way 
to directly see that the gas is superheated is to suddenly in- 
crease the coupling of the BEC to the thermal bath. In Fig. [3] 
we show the results of two experimental series in which a is 
quenched (within 10 ms) from 3 ao to 62 ao at different times 
in the superheated regime. The solid (open) symbols show 
iVo measured before (after) the quench. The small To is es- 
sentially unaffected by the change in a, and the sudden No 
decay is due to the increase in k [see Fig.[TJb)]. For reference, 
the green line shows the calculated N^ q at 3ao and orange 
shading indicates the superheated regime. 

As shown in Fig. [4j we have explored the limits of super- 
heating for a range of interaction strengths, including small 
negative values of a. For a < 0, a BEC is stable against col- 
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FIG. 4: Limits of superheating, (a) The highest temperature at which 
we observe a BEC, T, scaled to the equilibrium T c . Close to a — 
the BEC survives up to « 1.47T C . The red line shows the results of 
our numerical calculations, with the shaded area indicating the theo- 
retical uncertainty. Experimental error bars are statistical. The point 
at 62 ao is fixed to unity by the absolute atom number calibration 
1 191 1261 . (b) Temporal phase diagram. For each value of a we plot 
the equilibrium t c (green points) and the time tat which the BEC ac- 
tually vanishes (red points). The i errors correspond to dashed lines 
in Fig.[2]and the uncertainty in t c is indicated by the scatter of points. 
Solid curves are spline fits to the data. For a ~ the BEC survives 
in the superheated regime for a whole minute. Inset: numerically 
calculated phase diagram, with i data overlaid. 



lapse only for iVo < —C/a, with C « 2 x 10 4 ao for our 
trap parameters [23-25 1. However, after iVo drops below this 
critical value, at small \a\ it decays slowly. 

In Fig. |4ja) we plot the highest temperature at which we 
still observe a BEC, T = T(t — t), scaled to the equilibrium 
T c at the same N. For a —> 0, the BEC survives up to T w 
1.47 T c . (For comparison, this is analogous to superheated 
water at 275 °C.) 

In Fig. |4|b) we reconstruct the temporal phase diagram of 
our non-equilibrium gas. Here, a horizontal cut through the 
graph corresponds to a time series such as shown in Fig. [2] For 
each a, we plot the measured t (red points) and the equilibrium 
t c (green points). The solid curves are spline fits to the data. 
The width of the orange-shaded region corresponds to the time 
that the BEC survives in the superheated regime. For a w 
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this region spans a whole minute. 

The phase diagram in Fig.Qb) is measured by always start- 
ing with No ~ 2 x 10 4 . In general, non-equilibrium behaviour 
can strongly depend on the initial conditions. However, we 
find that t is essentially constant (within experimental errors) 
for initial Nq in the range (1 — 5) x 10 4 . The primary reason 
for this is that the three-body contribution to To grows with 
Nq; this leads to "self-stabilisation" of the condensed atom 
number on timescales much shorter than t. 

We theoretically reproduce our non-equilibrium observa- 
tions using a two-component model directly corresponding to 
Fig. |TJb). Starting with the measured initial No, we numer- 
ically simulate the evolution of a BEC coupled to a thermal 
bath characterised by To do this we calculate T from 

our experimental parameters, and for n we use the form G71 



(1) 



Here j e i oc a 2 is the elastic collision rate and A is a dimen- 
sionless coefficient. The largest uncertainty in our calcula- 
tions comes from the theoretical uncertainty in A « 1 — 10 
[28 1 . (For details see Supplementary Information 1 16 1.) 

In Fig. |4^a) we show the calculated T/T c . The red line 
corresponds to A = 3 and the shaded area to the range 
A = 1 — 10. The calculation generally captures our exper- 
imental observations well. With A = 3 we obtain quantitative 
agreement with the data, except exactly at a = 0, where the 
model is not valid. In the inset of Fig.|4|b) we show the calcu- 
lated temporal phase diagram, with A = 3, together with the 
experimental t data. Again the general features of the diagram 
are captured well for a ^ 0. 

In conclusion, we have observed superheating in a Bose- 
condensed gas with tuneable interactions, mapped out a non- 
equilibrium phase diagram of this system, and reproduced our 
measurements in numerical simulations based on a two-fluid 
picture of a partially condensed gas. The success of our cal- 
culations supports a conceptually simple way to think about 
dynamical non-equilibrium effects near a continuous phase 
transition. Extending the BEC lifetime by tuning interactions 
could also have practical benefits for precision measurements 
and quantum information processing. 



[1] Polkovnikov, A., Sengupta, K., Silva, A. & Vengalattore, M. 
Colloquium: Nonequilibrium dynamics of closed interacting 
quantum systems. Rev. Mod. Phys. 83, 863-883 (2011). 

[2] Kinoshita, T., Wenger, T. & Weiss, D. S. A quantum Newton's 
cradle. Nature 440, 900-903 (2006). 

[3] Winkler, K. et al. Repulsively bound atom pairs in an optical 
lattice. Nature 441, 853-856 (2006). 

[4] Sadler, L. E., Higbie, J. M., Leslie, S. R., Vengalattore, M. & 
Stamper-Kum, D. M. Spontaneous symmetry breaking in a 
quenched ferromagnetic spinor Bose-Einstein condensate. Na- 
ture 443, 312-315 (2006). 



[5 

[6 
[7 

[8 
[9 

[10 
[11 

[12 
[13 

[14 

[15 

[16 
[17 

[18 

[19 

[20 
[21 

[22 

[23 

[24 

[25 
[26 

[27; 
[28 



Hofferberth, S., Lesanovsky, I., Fischer, B., Schumm, T. & 
Schmiedmayer, J. Non-equilibrium coherence dynamics in one- 
dimensional Bose gases. Nature 449, 324-327 (2007). 
Haller, E. et al. Realization of an excited, sU'ongly correlated 
quantum gas phase. Science 325, 1224-1227 (2009). 
Guzman, J. et al. Long-time-scale dynamics of spin textures in 
a degenerate F = 1 87 Rb spinor Bose gas. Phys. Rev. A 84, 
063625 (2011). 

Cheneau, M. et al. Light-cone-like spreading of correlations in 
a quantum many-body system. Nature 481, 484-487 (2012). 
Trotzky, S. et al. Probing the relaxation towards equilibrium in 
an isolated strongly correlated one-dimensional Bose gas. Na- 
ture Physics 8, 325-330 (2012). 

Mark, M. J. et al. Preparation and spectroscopy of a metastable 
Mott-insulator state with attractive interactions. Phys. Rev. Lett. 
108,215302 (2012). 

Smith, R. P., Beattie, S., Moulder, S., Campbell, R. L. D. & 
Hadzibabic, Z. Condensation dynamics in a quantum-quenched 
Bose gas. Phys. Rev. Lett. 109, 105301 (2012). 
Gring, M. et al. Relaxation and pre-thermalization in an iso- 
lated quantum system. Science 337, 1318-1322 (2012). 
Mathey, L. & Polkovnikov, A. Light cone dynamics and re- 
verse Kibble-Zurek mechanism in two-dimensional superfluids 
following a quantum quench. Phys. Rev A 81, 033605 (2010). 
Fedichev, P. O., Shlyapnikov, G. V. & Walraven, J. T. M. Damp- 
ing of low-energy excitations of a trapped Bose-Einstein con- 
densate at finite temperatures. Phys. Rev. Lett. 80, 2269 (1998). 
Pethick, C. & Smith, H. Bose-Einstein Condensation in Dilute 
Gases (Cambridge University Press, Cambridge, 2002). 
Online supplementary information. 

Campbell, R. L. D. et al. Efficient production of large 39 K 

Bose-Einstein condensates. Phys. Rev. A 82, 063611 (2010). 

Roati, G. et al. 39 K Bose-Einstein condensate with tunable 

interactions. Phys. Rev. Lett. 99, 010403 (2007). 

Smith, R. P., Campbell, R. L. D., Tammuz, N. & Hadzibabic, Z. 

Effects of interactions on the critical temperature of a trapped 

Bose gas. Phys. Rev. Lett. 106, 250403 (2011). 

Tammuz, N. et al. Can a Bose gas be saturated? Phys. Rev. 

Lett. 106, 230401 (2011). 

Smith, R. P. & Hadzibabic, Z. Effects of interactions on 
Bose-Einstein condensation of an atomic gas. arXiv: 1203.2063 
(2012). 

Dalfovo, F, Giorgini, S., Pitaevskii, L. P. & Stringari, S. The- 
ory of Bose-Einstein condensation in trapped gases. Rev. Mod. 
Phys. 71, 463-512 (1999). 

Ruprecht, P. A., Holland, M. L, Burnett, K. & Edwards, M. 
Time-dependent solution of the nonlinear Schrodinger equation 
for Bose-condensed trapped neutral atoms. Phys. Rev. A 51, 
4704^1711 (1995). 

Gerton, J. M., Strekalov, D., Prodan, I. & Hulet, R. G. Direct 
observation of growth and collapse of a Bose-Einstein conden- 
sate with attractive interactions. Nature 408, 692-695 (2000). 
Donley, E. A. et al. Dynamics of collapsing and exploding 
Bose-Einstein condensates. Nature 412, 295-299 (2001). 
Smith, R. P., Tammuz, N., Campbell, R. L. D., Holzmann, M. 
& Hadzibabic, Z. Condensed fraction of an atomic Bose gas 
induced by critical correlations. Phys. Rev. Lett. 107, 190403 
(2011). 

Gardiner, C. W, Zoller, P., Ballagh, R. J. & Davis, M. J. Ki- 
netics of Bose-Einstein condensation in a trap. Phys. Rev. Lett. 
79, 1793-1796 (1997). 

Gardiner, C. W., Lee, M. D., Ballagh, R. J., Davis, M. J. & 
Zoller, P. Quantum kinetic theory of condensate growth: com- 
parison of experiment and theory. Phys. Rev. Lett. 81, 5266- 



5 



5269 (1998). 

* These two authors contributed equally. 
t rps24@cam.ac.uk 
Acknowledgements 

We thank S. Beattie and S. Moulder for experimental as- 
sistance. This work was supported by EPSRC (Grant 



No. EP/K003615/1), the Royal Society, AFOSR, ARO and 
DARPA OLE. 
Author Contributions 

All authors contributed extensively to this work. 
Additional Information 

The authors declare no competing interests. Correspondence 
and requests for materials should be addressed to R. P. S. 



Supplementary Information 



6 



Determination of TV, E and No 

We take an absorption image of the atom cloud after r = 18 ms of time-of -flight (TOF) expansion from a nearly isotropic trap 
with a geometric mean of trapping frequencies u>/2ir ~ 70 Hz. 

For absolute calibration of our atom numbers we use a T c measurement at a = 62 ao, assuming that at this a the cloud is in 
equilibrium (HE). 

For a < 100 ao we assess the interaction-energy contribution to the total energy E to be < 1% and thus E w 2E k , where 
Ej. is the kinetic energy. We obtain Ek from the second moment of the atom distribution after TOF, and correct it for the small 
effect of the initial in-trap cloud size. This amounts to rescaling the energy by a factor w 2 r 2 /(l + uj 2 t 2 ). 

In Fig. S|5]we show TV and Ef. for the same two experimental series shown in Fig. 2 in the main text. 

To improve the detection of small TVq values we always switch a to zero at the start of TOF [ 1 1. 
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Fig . S5 : Parametrization of TV and Ek ■ The total atom number TV and kinetic energy per particle Ek /TV for two experimental series at 83 ao 
and 5 ao . The experimental data are fitted with polynomial forms to obtain smooth functions TV (t ) and Ek (t) . 



Equilibrium calculations 

For an equilibrium ideal Bose gas in a spherically symmetric harmonic trap of frequency cj, the thermal component satisfies: 



E k (p',T) = 3 /o°°g2 (cxp(4j#-^)) jx&z 

where TV° = ((3) (^|j-) 3 ^1 — §$5) ~k^T ) includes the finite size correction, — ^Tluj, and g2(z) is the dilogarithm 

function. For TV > TV,?, the chemical potential is capped at // = /i° and also fiQ = The thermal atom number is saturated, 
TV' = TV?, and any additional particles must go into the condensate, TV = TV — iV°. 
Interactions modify this picture in two ways: 



7 



(i) The critical point is shifted. At mean-field (MF) level \i c = fx® + 4£(3/2)a/A, where A = J is the thermal 
wavelength. A small beyond-MF correction is quadratic in a/ A and has an additional logarithmic correction. Experimentally, 
the corresponding N c is (HE) 

N c = 7Y C ° ^1 - 3.426^ + 42 (|)^ . (4) 

(ii) Due to interactions, in the presence of a BEC, N' is no longer saturated at N c . Empirically, 

N> = N c + So(N ) 2 / 5 + S 2 (N ) i/5 , (5) 

where the non-saturation parameters Sq and S% depend on a and T J3l S) . The excess number of thermal atoms, N' — N c , can be 
directly attributed to the shift of the chemical potential above /i c ; for an interacting BEC [1q > fi c and in equilibrium // = /j,q. 
For jUo we use a modified Thomas-Fermi law |5 1: 



1 '15^ ■ r -'V"'"' 3\. ,6, 

O'osc 



where a osc = ^Jh/mui is the harmonic oscillator length. 

Eq. (|6ji and Eqs. Q and ^ modified to include interaction effects (N° — > N c and — > f^c) form a complete set needed for 
our calculations |6|. We proceed in two ways: 

(1) Assuming that the system is in global equilibrium, fjf = fio = fj, cq , we use only N(t) and E(t) to calculate /J. eq (t), NQ q (t) 
(green lines in main-text Fig. 2) and T cq (t). 

(2) Additionally, from N, E and the measured Nq we calculate (iq and fi' without assuming // = /j,q (red and blue points in 
the bottom panel of main-text Fig. 2). 



Justification of the two-fluid picture for a gas out of global equilibrium 

In our theoretical picture (Fig. 1(b) in the paper), we assume that in the superheated regime the thermal and condensed 
components can still to a good approximation be assumed to be separately in equilibrium. Moreover, we assume that they are at 
the same temperature, but just have different chemical potentials. 

Here we provide a more detailed justification of these assumptions. 

First, for the thermal component we show in Fig. S|6]that the radial velocity distribution still looks like a thermal distribution 
at a temperature very close to the calculated T cq (N, E). Here we show data for the same 5 ao series as shown in Fig. 2(b) in 
the main paper. In Fig. S|6|a) we show the distribution measured at t = 45 s, i.e. deep in the superheated regime. The data 
(red) is fitted almost perfectly by an equilibrium thermal distribution constrained to be characterised by N, E and T cq (green). 
An unconstrained fit (blue) gives only a very slightly different shape with within few % of T cq . In Fig. sj^Jb) we compare 
the calculated T cq (green) and the fitted T* (blue) for the whole 5 ao series. For comparison we also show the equilibrium T c 
(black line). Note that this is the same plot as in Fig. 2(b) in the main paper, with just the T* points added. 

Second, for the collective excitations in the BEC to be in equilibrium with the thermal bath, the Landau-damping time tl [7 8] 
must be short compared to the characteristic time scale of the experiment. For a uniform system at a temperature higher than the 
interaction energy per particle |[8l : 



1/2 -v 2 

tlu _ V A^ 

2tt ~ 4 7 ra 1 / 2 ' ^ ' 

where uj is the excitation frequency, uq is the condensate density, and A^ is the thermal wavelength. For our assessment we use 
this uniform-system result with our peak no; this only overestimates tl for a harmonically trapped gas Q. 

Exactly at a — the damping time diverges and our theoretical picture breaks down. However, already for a — 1 a , for all 
our experimental parameters the RHS of Eq. (j7| is < 100. Then, even for our lowest-energy modes, with lo/2tt ~ 100 Hz, we 
get t l < 1 s. 
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Fig. S6: Thermal distribution in a gas out of global phase equilibrium, (a) For a gas at 5 ao in the superheated regime we show the radial 
velocity distribution (red), the distribution corresponding to N, E and T cq (green) and the unconstrained fit (blue) giving T. Even though the 
gas is not in true global equilibrium the distribution still looks thermal and T' and T oq agree to within few %. (b) Comparison of T s (blue) 
and T cq (green) for the whole 5 ao data series shown in Fig. 2(b) of the main paper. The solid black line shows the equilibrium T c . 



Non-equilibrium evolution of No 



The non-equilibrium evolution of No is described by the differential equation 

Nq = -k- T N , 



(8) 



where k is the coupling to the thermal bath due to elastic collisions and Tq is the instantaneous inelastic loss rate per particle. 
Following |9| we use 



K = A 



8m(ak B T) 2 2j 3(//-^) 



,/3(A"o-/*') 



1 



N = Aj el N 



(9) 



__ 8m(ak B T) 2 



is essentially the elastic scattering rate for a thermal cloud at \i! and Am 1 — 10 is a 



where 7cl - ^ 
theoretically uncertain pref actor Q. 

In the inelastic loss term we include contributions from one-body scattering and three-body recombination, Tq — + T 
In our system, one-body loss is dominated by spontaneous scattering of photons from the trapping laser beams. We calculate it 
from the known wavelength and intensity of the beams. For all the reported experiments 
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The loss rate of condensate atoms due to three-body recombination in the presence of a thermal cloud is given by ifTUl 

r( 3 > = ^(< n 3)+6<non0 + 6<n' 2 >) , 
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where n is the condensate density, n' the thermal density, K 3 (a) the known a-dependant three-body coefficient ATI [12 1. and 
(...) stands for an average over the density distribution. We set n' to its value in the centre of the trap and for the condensate we 
again use a modified Thomas-Fermi approach: 
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) 3 / 2 . Eqs. ( 12 1 and ( 13 1 smoothly interpolate between the ground state result (for 

1). Note that other forms which smoothly interpolate 



Here GS refers to the non-interacting Gaussian ground state and TF to the Thomas-Fermi approximation, (uo)tf 
V2^(15/7)(15AWaoscr 3/5 (no) G s and (n ) TF ^ V3Tril5 2 /56)(15N a/a osc )- 6 / 5 {n$) GS , where (n ) GS 
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N a/a osc -C 1) and the Thomas-Fermi approximation (for N a/a, 
between these two limits give essentially the same results. 

Using Eqs. ^ - ( 13 i we simulate the evolution of the condensate atom number, No(t), from its initial value AT (i = 0). We 
use the measured N(t) and E(i) and the numerically evolved No(t) to obtain fi'(t) for use in Eq. &j. To determine t from our 
calculations we define the condensate to be present if No is larger than A*"™ 1 " = 3/cb? 1 / (hw), the thermal occupation of the first 
excited state. 

In the main text we show the results of our calculations for T and t (Fig. 4), with A — 3. In Fig. SjT] we show that our 
calculations (with the same value of A) also describe well the full dynamics No(t). 
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Fig. S7: Non-equilibrium No dynamics. We plot the calculated No(t) (dashed red line) together with the measured iVo (red points) for the 
same 5 ao data series as in Fig. 2(b) in the main text. For comparison we also show the calculated iVg q (f) (solid green line). 
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